An explicit finite-difference scheme for simulation of moving particles

We present an explicit finite-difference scheme for direct simulation of the motion of solid particles in a fluid. The method is based on a second order MacCormack finite-difference solver for the flow, and Newton’s equations for the particles. The fluid is modeled with fully compressible mass and momentum balances; the technique is intended to be used at moderate particle Reynolds number. Several examples are shown, including a single stationary circular particle in a uniform flow between two moving walls, a particle dropped in a stationary fluid at particle Reynolds number of 20, the drafting, kissing, and tumbling of two particles, and 100 particles falling in a closed box. Comments Postprint version. Published in Journal of Computational Physics, Volume 212, Issue 1, 2006, pages 166-187. Publisher URL: http://dx.doi.org/10.1016/j.jcp.2005.06.021 This journal article is available at ScholarlyCommons: http://repository.upenn.edu/meam_papers/72 An Explicit Finite-Difference Scheme for Simulation of Moving Particles A. Perrin, H. H. Hu Department of Mechanical Engineering and Applied Mechanics University of Pennyslvania 297 Towne Building, 220 S. 33rd Street, Philadelphia, PA 19104 Abstract We present an explicit finite-difference scheme for direct simulation of the motion of solid particles in a fluid. The method is based on a second order MacCormack finitedifference solver for the flow, and Newton’s equations for the particles. The fluid is modeled with fully compressible mass and momentum balances; the technique is intended to be used at moderate particle Reynolds number. Several examples are shown, including a single stationary circular particle in a uniform flow between two moving walls, a particle dropped in a stationary fluid at particle Reynolds number of 20, the drafting, kissing, and tumbling of two particles, and 100 particles falling in a closed box.We present an explicit finite-difference scheme for direct simulation of the motion of solid particles in a fluid. The method is based on a second order MacCormack finitedifference solver for the flow, and Newton’s equations for the particles. The fluid is modeled with fully compressible mass and momentum balances; the technique is intended to be used at moderate particle Reynolds number. Several examples are shown, including a single stationary circular particle in a uniform flow between two moving walls, a particle dropped in a stationary fluid at particle Reynolds number of 20, the drafting, kissing, and tumbling of two particles, and 100 particles falling in a closed box.